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We investigate the vortex structure in chiral p-wave super- 
conductors by the Bogoliubov-de Gennes theory on a tight- 
binding model. We calculate the spatial structure of the 
pair potential and electronic state around a vortex, includ- 
ing the anisotropy of the Fermi surface and superconducting 
gap structure. The differences of the vortex structure between 
sinpx+isinpji-wave and smp x — isinp^-wave superconductors 
are clarified in the vortex lattice state. We also discuss the 
winding ^3 case of the sin (p x + p y ) ± isin(— p x + p y )-wave 
superconductivity. 

PACS numbers: 74.60.Ec, 74.20.Rp, 74.70.Pq, 74.25.Jb 



I. INTRODUCTION 

Recently, much attention has been focused on 
the superconductivity in quasi-two-dimensional metal 
Sr2Ru04,EJ there is experimental evidence that it anight 
be a realization of a chiral p-wave superconductor .0 The 
17 0-NMR measurements reported that there is no re- 
duction of the Knight shift in the superconducting state, 
supporting the spin triplet pairing!! Moreover /iSR mea- 
surements claim that spontaneous magnetic moment ap- 
pear in the superconducting state, suggesting a pairing 
state with broken time reversal symmetry.cl Therefore, 
the symmetry of superconductivity is likely to be that of 
the chiral p-wave with the basic form A± ~ p x ± vp y and 
inplane equal-spin pairing. This simplified form shall, 
however, not preclude any details of the momentum de- 
pendence of the quasiparticle gap. Detailed discussions 
on the pairing symmetry were given in Refs. H [|. 

At zero field, A + and A_ state are degenerate and, 
in general, domain formation of the two states, A + and 
A_, is expected. This degeneracy is removed under an 
external magnetic field with a component perpendicular 
to the basal plane, since A± corresponds to states with 
an orbital angular momentum along the z-axis. Also the 
vortex structure in the mixed phase is different for A + 
and A_ and we may consider the so-called positive vor- 
tex (P-vortex) and negative vortex (N- vortex) £3 In the 
former (latter) corresponds to a vortex with winding ori- 
entation parallel (antiparallel) to the internal winding of 
the Cooper pair. If the external field H is parallel to z 
for e > (2e is the charge of the Cooper pair) or, if H is 
antiparallel to z for e < 0, the vortex with a winding +1 
appears, i.e., A(r) ~ /(r)e 1<6 where r = (r cos <fi, r sin <fi) . 



This corresponds to the P-vortex (N- vortex) for the phase 
A + (A_). For reversed magnetic field direction, one finds 
vortices with opposite winding —I, thus A(r) ~ /(r)e _Ic ^, 
which represents the N- vortex (P-vortex) for A + (A_). 

The differences of the structure between P-vortex and 
N-vortey. were studied within the GinzburgpLa i idau (GL) 
theoryOO and the quasiclassical theory.EijO Further- 
more the electronic state around the vortex core was 
studied by the Bogoliubov-de Gennes (BdG) Xheopji in 
the single vortex case and continuum model. E2rli3li3 It 
was suggested that the character of the quantized energy 
level for bound state of quasiparticles around the vortex 
core is different for P-vortex and N-vortex. 

In this paper we investigate the difference of the vor- 
tex structure between the P-vortex and N-vortex, based 
on the BdG theory for a tight-binding model, consider- 
ing the form of pair potential and electronic states in the 
vortex lattice state. So far, this method has been used 
to the study of the c4,2_ y 2-wave superconductivity in the 
high-T c su]?€jrcanductors, applied— tp-ithe extended Hub- 
bard modelEjila or the t-J model.Ejc3 Here, we introduce 
the spin triplet p-wave pairing interaction instead of the 
singlet d-wave and s-wave pairings. 

Using the tight-binding model allows us to take cer- 
tain aspects of the real band structure into account. 
While Sr2Ru04 is a metal with three electronic bands, 
we will restrict to one band only. This is justified from 
the point of view that one of the three bands domi- 
nates the superconducting propertiesBEJ We take the 
nearly cylindrical symmetric band originating fronp-the, 
Ru 4+ -4d t2g orbital d xy yielding the so-called 7-band.E3~0 
The square lattice tight-binding model includes nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) hop- 
pingjUsing the parametrization given, for example, in 
Ref. Eq. Any dispersion along the z-axis is neglected. The 
pairing interaction is taken on the lattice as real space 
attractive NN interaction yielding the gap function form 
sinp x ± i sinpj,£3'c3 The superconducting gap is nodeless 
on the Fermi surface, but has a strong anisotropy. 

Our self-consistent calculation of the vortex structure 
reveals the anisotropy of the gap and the band struc- 
ture, when we consider the local density of states (LDOS) 
around the vortex core. Moreover, we observe the de- 
tailed structure of the pair potential in the vicinity of 
the vortex core, which contains induced components of 
A_ (A + ) in the background of dominant A + (A_) super- 
conductivity. Since our discussion includes the complete 
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vortex lattice, we are able to discuss the contribution 
of the inter-vortex transfer of low-energy quasiparticles 
bound at the vortex core. 
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FIG. 1. (a) Two dimensional Fermi surface used in our cal- 
culation. We set parameters as t' = 0A7t and /i = 1.2i, which 
are appropriate to the 7-sheet in S^RuCU. In the phase of the 
superconducting gap function A+(p) = shipa, + isinpy and 
A' + (p) = sin^ +p y ) + isin(-p a; +p y ), the F point gives a 
winding +1. Futher, A+(p) has four additional singularities 
at (±7r/2, ±7r/2)(closed circles in the figure), which give each 
a winding — 1. (b) Fermi velocity |vF(p)|(arbitrary unit) and 
the amplitude of the gap functions |A+(p)| (solid line) and 
|A' + (p) I (dotted line) along the Fermi surface line. 9t F = 
(= 7r/4) corresponds to T-M (r~X) direction, (c) The di- 
rectional dependent coherence length £(k) = |ir(p)|/|A(p) 
along the Fermi surface line. Solid (dotted) line is for the 
|A+(p)| (|A' + (p)| ) case, (d) The phase of A+(p) along the 
Fermi surface line. We also show ReA+(p) and ImA+(p). (e) 
The phase of A' + (p) along the Fermi surface line. 



Furthermore we consider the case of NNN pairing inter- 
action, in contrast to NN pairing mentioned above. The 
gap function then corresponds to the form sin (p x +p y )± 
i sin (— p x +p y ). We will see that this case gives =f3 wind- 
ing of the Cooper pair on our considering Fermi surface 
and belongs therefore to a different topological class. In 
our study we will see the role of the Cooper pair winding. 

After describing our formulation of the BdG theory 
on tight-binding model in Sec. ||, we explain the pair 
potential structure on the Fermi surface in Sec. |!Tl[ Wc 
consider the pair potential structure of the vortex lattice 
state in Sec. |Fv|, and the electronic state in Sec. [v|, 



comparing P-vortex and N-vortex both in the NN site 
pairing case and the NNN site pairing case. The last 
section is devoted to the summary and discussion of our 
results. 



II. BOGOLIUBOV-DE GENNES THEORY ON 
TIGHT-BINDING MODEL 

We begin with the tight-binding model on a two- 
dimensional square lattice of Ru-site, and introduce a 
pairing interaction between NN sites or NNN sites. The 
Hamiltonian for the system in a magnetic field is given 
by 



pair 



(2.1) 



with the creation (an annihilation) operator a\ (<Zj, a ) for 
spin a =f, I and site i = (i x ,i v ). The transfer integral 
is expressed as 



= U.j exp[i— / A(r) • dr] 



(2.2) 



with the vector potential A(r) = x r in the symmet- 
ric gauge, and the flux quantum <po. For NN-hopping, 
tij = t and for NNN-hopping in diagonal direction of 
the square lattice plaquette, t% j — t' . To reproduce the 
Fermi surfacer-topology of the 7-sheet in Sr2Ru04, we 
set t' — 0.47i,c3 in order to obtain a closed electron-like 
Fermi surface for a filling of n ~ 1.12, i.e., beyond half 
filling. Keeping this condition, our results do not vary 
qualitatively with different parameters. 

The spin part of pairing interaction between i- and j- 
sites is decomposed to the singlet component go.ji and 
triplet components g x ,ji, 9y,jii 9z,ji as follows, 



Impair — 2 



i,j,ai~a4 m— 0,x : y : z 



with the Pauli matrices a x , a y , a z and a unit matrix uq. 
The subscripts a\ ~ are spin indices. If we consider 
the following interaction 
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^pair 2 ^ ^ (Vji'fl'jTli Jx.jiSx J Sx ,i 

~^~Jydi^y>3^'y>i ^z,ji^z,j^z,i) (^*^) 

using the number density operator nt and the spin den- 
sity operator S x> i, S y ^, S z ,i at each site, the pairing in- 
teractions are obtained as 

%90,ji — Vji — (Jx,ji + Jyji + Jz,ji)/^i (2-5) 

*^9x,ji Vji ~t~ ( Jx,ji ~t~ Jy,ji J z.ji) / (2-6) 
^9y,ji Vji i^x.ji Jy,ji ^Z,ji) (2-7) 

^9zji — ^ji ~t~ {Jxji ~t~ Jyji d z ji)/4i. (2-8) 



The superconducting order parameter is decomposed 



to 



(2.9) 



The self-consistent condition for the ci-vector is written 
as 



d m ,ji — 9m,ji ^ ] ( (J mi< J y)a 3 a 1 ( a ia 1 a ja 3 ) ■ 



(2.10) 



In this paper, we set goji = g x ,ji = 9y,ji — so that 
do.ji — d x ,ji — dy.ji — 0, since we consider the case of the 
spin triplet pairing with d z only. Thus, when we assume 
the pairing interaction works only between NN (NNN) 
sites, we set g z ji — g z for NN (NNN) site pairs, and 
otherwise g z ji = 0. 

In terms of the eigen-energy E e and the wave func- 
tions u e (ri), v e (ri) at i-site, the Bogoliubov-de Gennes 
equation is given by 
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K d ■■ 



where K a ^ 
eigenstates.ll 



u e (ri) 



= E f 



u e (r 3 ) 

Ve{Vj) 



, (2.11) 



/J,5ij and e is an index of the 



The self-consistent condition in Eq. ( 2. 1C ) is reduced 



to 



with 



d z ,ji = 9z.ji((ajia. n ) + (a jT a a )) (2.12) 

(a jiai1 ) = '£v* e (r j )u £ (r i )f(E e ), (2.13) 

e 

(a n a tl ) ='£u € (r j X(r i )f(-E e ), (2.14) 



and the Fermi distribution function f(E). The spin- 
triplet pair potential satisfies the relation d z ^j — —d z ji, 
i.e. it corresponds to odd-parity pairing. 

The orbital part of the pair potential at each site i can 
be decomposed into a sinp x and a sinp y component as 



Ap x (r,) = (A*,; - A_ 4li )/2, 
A Pa (r 4 ) = (A Si< - A_ fi|i )/2 



(2.15) 
(2.16) 



in the NN site pairing case, or into a sin (j> x + p y ) and a 
sin (— p x + p y ) component as 



^Px+Py( r i) — i^x+y.i A_£_y j i)/2, 
Px~\~Py\^i) — &-\-yA Ax— y,i)/2 



(2.17) 
(2.18) 



in the NNN site pairing case. In Eqs. ( |2.15| )-( |2.18| ) we 
denoted 



Ag,i = d Z:iii+ e exp[i- 



(n+r i+s )/2 



A(r)-dr]. (2.19) 



When the pair potential is uniform, our BdG formu- 
lation is reduced to the conventional theory for p- 
wave superconductors with the pairing functions sinp x 
and sin p y , or sin (p x + p y ) and sin (— p x + p y ), by the 
Fourier transformation to the momentum space. For 
sin Pj, ± isinpy-wave superconductivity, we define the 
pair potential as A±(r^) = A P;r (rj) =F iA PH (r i ), and 
for sin (p x + p y ) ± i sin (— p x + p y )-w&ve superconductiv- 
ity, A' ± (ri) = A Pa+Py (ri) t iA_ Pa , +P!/ (ri). 

To investigate the electronic structure around the vor- 
tex, we calculate the LDOS N(E,n) = N^(E,n) + 
Ni(E,ri) at the i-site, where 

N^E,n) = \uM)\ 2 S(E - E € ), (2.20) 
AT;( J B,r i ) = X)KN| 2 ^ + J Be) (2.21) 

e 

for up-spin and down-spin contributions, respectively. 
The electron number density at each site is given by 
n(r-j) = n T (ri) + n i (r l ) with 

n T (r,)= f dEN^E,n)f(E) =J2We(^)\ 2 f(Ee), (2.22) 

n l (r i )= [ AEN^E^fiE) = Y J W{v l )?f{-E i ). (2.23) 

We consider a system with a square unit cell of N r x N r 
sites, where two vortices are accommodated. The NN 
vortices of the square vortex lattice are located at the 
45° directions from a-axis of the atomic structure. This 
vortex lattice configuration-js suggested from the neu- 
tron scattering experiment .E3 By introducing the quasi- 
momcntum of the magnetic Bloch state 



u e (r) = u k: e k (r)e 1 
with 

2?r 



v e (r) = £ k ,e k (r)e 



iki 



(2.24) 



cN r N k 



(l x ,l y = l,---,N k ) (2.25) 



we obtain the wave function under the periodic bound- 
ary condition whose region covers Nj~ x Nk unit cells (c 
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is a lattice constant). The eigenstate e of Eq. (2.11) 
can then be labelled by the q uasi- momentum k and the 
eigenenergy ek based on Eq. ( 2.24 ). 

Our calculation starts from the initial state of a square 
vortex lattice solution v I / o( r ) m the lowest Landau level. 
Thus, when the dominant superconductivity is A± in 
the NN-site pairing case, the initial state is given by 
(A fe fa),A p „fa)) = (l,±i)*o(r). A' ± in the NNN 
site pairing corresponds to (A p „ +P „ fa;), A_ p „ +P „ fa)) = 
(1, ±i) 1 I f o(r). We iterate the calculation of Eqs. (2.11)- 
(2.14) until convergence is achieved. We typically con- 
sider the case N r = 20 and g z ji = — l.Ot. The spatially 
averaged electron density is set to n = nfa) ~ 1.12 by 
tuning the chemical potential ji. We consider a sequence 
of temperatures between T — and T c . With our choice 
of parameters, we find T c ~ 0.27i. In the figures of this 
paper, we mainly show the case for T = O.lt. We con- 
sider two cases for the position of the vortex center. In 
one case, the vortex center is located just on the atomic 
site. We call it as site-centered vortex. For the other case, 
it is located in the middle of square plaquette surrounded 
by four atomic sites. We call it as plaquette-centered vor- 
tex. We mainly focus the site-centered vortex case. We 
obtain qualitatively the same pair potential structure and 
the electron number density structure in both cases. 



III. PAIR POTENTIAL STRUCTURE ON THE 
FERMI SURFACE 

Before discussing the vortex structure, we analyze the 
shape of the Fermi velocity and the superconducting 
gap function on the Fermi surface at zero field. The 
two dimensional Fermi surface of our model is shown 
in Fig. 0(a), corresponding to the 7-sheet of Sr2Ru04. 
Along this Fermi surface line, the angle dependence of 
the Fermi velocity v-p — |vp| is depicted in Fig. 0(b). 
It has a maximum for 9k F = (T-M direction) and 
minimum at 6k F — tt/A (T-X direction), which results 
from the vicinity of a van Hove singularity in the T- 
M direction. In the same figure also the magnitude 
of the gap function A + (p) = smp x + ismp y on the 
Fermi surface is plotted. |A+(p)| has minimum at 
0k F — 0. The resulting directional dependent coher- 
ence length £(6> fcF ) = t> P (0 fep )/|A + (0 fep )| is shown in 
Fig. 0(c), where £ is largest for 9k F — 0, and the ratio 
£(#fc F = 0)/£(#fc F = 7r/4) is nearly 5. The directional de- 
pendence of £ (8k F ) is important for the qualitative inter- 
pretation of the spatial distribution of low energy quasi- 
particles around the vortex, as discussed in Sec. [v|. For 
A+(p) = sin fac + p y ) + isin(— p x +p y ), the amplitude 
|A^_(p)| has minimum at 9k F — 0. While it has local 
minimum at 9k F — 7r/4, it is larger than \A' + (9k F = 0)|. 
Therefore, £(#fc F ) has maximum at 9k F = also for this 
A' + (p). The ratio £(6» feF = O)/£(0 fep = tt/4) is decreased 
to 2.7. 

In the chiral p-wave superconductor, it is important 



to consider the phase winding of the gap function on 
the Fermi surface. In Fig. 0(d), we present the phase 
argA + (#fc F ) with ReA + (6'fc F ) and lmA + (9k F )- The phase 
of A + (9k F ) undergoes the winding of +1 (x2tt), when we 
follow it around the Fermi surface. On the other hand, 
the phase of A' + (9k F ), as shown in Fig. 0(e), acquires 
a winding —3 around the Fermi surface. The origin of 
this difference lies in the structure of the singularities 
in the two gap functions. The gap function A+(p) and 
A' + (p) vanish at the T point, where their phase shows 
a +1 winding. Furthermore, A^_(p) has additional zeros 
at the four points (p x ,p y ) = (±7r/2, ±7r/2) not present in 
A + (p) . These additional singularities of the gap function 
give each a winding — 1. The Fermi surface in Fig. 0(a) 
is enclosing all five singularities of A' + (p) (a +1 winding 
point at r and four —1 winding points at (±7r/2, ±7r/2)). 
Hence, the total winding of A' + (p) is —3. A sufficiently 
small Fermi surface centered around the T-point, that 
does not enclose these outer singularities would, conse- 
quently, have only a winding +1, like A + (p). This differ- 
ence of the winding structure around the Fermi surface 
affects the low energy quasiparticle structure around the 
vortex core. 



IV. PAIR POTENTIAL STRUCTURE IN THE 
VORTEX LATTICE STATE 

First, we consider the spatial structure of the domi- 
nant pair potential in the vortex lattice state. Figures 
0(a) and 0(b) show the amplitude of A+ for P-vortex 
and A_ for N-vortex in the case of vortex with winding 
+1, respectively. Naturally, the amplitude is suppressed 
around the vortex core. As seen from contour lines in 
Fig. 0, the vortex core has square-like shape, reflecting 
the fourfold symmetry of the Fermi surface structure and 
the superconducting energy gap | sinp^ ± isinpjj in the 
momentum space. At low temperature, the orientation of 
the square shape is different depending the winding. In 
the P-vortex (N-vortex) case, the corner of the square is 
directed to the 0°- (45°-) direction from the a-axis, since 
the amplitude is more suppressed in this direction. This 
difference is smeared with increasing temperature. At 
higher temperature near T c , the vortex core radius is in- 
creased, and the contour lines both for P- and N-vortices 
are reduced to the same structure as the Abrikosov vor- 
tex solution |\&o( r )|j where the amplitude is suppressed 
along the NN vortex direction. 
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FIG. 2. Amplitude of the dominant component of the pair 
potential for P-vortex (a) and for N- vortex (b) at T = O.lt 
for sin pa; ± isinp^-type superconductivity. We, respectively, 
plot |A + (r)|/|# z | and \A-(r)\/\g z \ in the area of 20 x 20 
sites, where vortices are located in the middle and at the 
four corners of the figure, (c) and (d), respectively, show 
the dominant component |AY(r)|/|<? z | and |A'_(r)|/|<? z | in the 
sin (p x +p y )±i sin (— p x + pjj-wave case. 




FIG. 3. Amplitude of the induced component of the 
pair potential for P-vortex (a) and for N-vortex (b) at 
T = O.lt for the sinpa: ± isinp^-type state. We, respec- 
tively, plot |A_(r)|/|<; z | and |A + (r)|/|g z |. The plotted area 
is the same as in Fig. |2| (c) and (d), respectively, show 
the induced component |A'_(r)|/|<; z | and |A+(r)|/|<7 z | for the 
sin (p x + p y ) ± i sin (— p x + p y )-type state. 

Next, we consider the additionally induced component 
of the pair potential. In the vortex core region where the 
dominant component of Fig. ^| shows strongest spatial 
dependence, the other pairing component is appears. It 
is known that this component occurs also in the frame- 
work of the GL theory, where A + and A_ couple to each, 
other through the gradient terms in the GL equation.Eil 
Figure || shows the amplitude of A_- (A+-) component 



in the dominant A + (A_) pairing channel. The induced 
component has the fourfold symmetric structure. This 
shape is different between P- and N-vortices. In the P- 
vortex case, the induced component |A_| is sharply sup- 
pressed along the line of 45° direction from a-axis. This 
difference comes from the winding structure of the in- 
duced component. At higher temperature near T c , this 
difference is smeared. 

The winding structure is schematically shown in Fig 
[|. For the P-vortex, Fig. ||(a) shows the basic feature of 
winding at low temperature. Fig. ^(b) shows the wind- 
ing behavior at higher temperature for the P-vortex and 
for any temperature for the N-vortex. First, we consider 
the N-vortex. The induced A + -component of Fig. ||(b) 
has winding —1 at the vortex center, where the dominant 
A --component has winding +1. Since the total of wind- 
ing number should be +1 per vortex both for the A + - and 
the A_-component, the A_|_-component has also winding 
+2 at the mid points between NNN vortices (i.e. each 
corner of the boundary line of the square vortex lattice 
unit cell). 

Now we consider the P-vortex case at low tempera- 
ture. There, the winding structure depends on temper- 
ature and applied magnetic field. If vortex core radius 
r coro is short enough compared with the inter-vortex dis- 
tance It, — ((fio/H) 1 / 2 at low temperature, the induced 
A --component has winding —1 at the vortex center, and 
winding —2 at the corner of the boundary line. Fur- 
thermore, there are four points with winding +1 in the 
45° direction near the vortex core. Since |A_| = at 
these winding points, the amplitude |A_| is sharply sup- 
pressed along the 45° direction, as shown in Fig. ||(a). It 
is known that the induced s-wave component in d x 2_ y 2- 
superconductor shaajSj-the same type suppression due to 
the winding points In the limit of an isolated single 
vortex (l v /r COIC — * oo), four +l-winding points approach 
the — 1-winding point at the vortex center. They turn 
together into a +3-winding point at the vortex center. 
This structure of winding +3 for P-vortex wa s obtai ned 
in the calculation for an isolated sing le vortex.BBt3 On 
the other hand, when r COTC /l v is increased with raising 
temperature or applied field, each winding +1 point ap- 
proaches the boundary line, i.e., the midpoint between 
NN vortices. After approaching the boundary, these 
winding +1 points are combined with winding —2 points 
at the corners of the boundary line and they become 
a +2-winding point at the corner. The combination of 
winding points easily occurs at the boundary, since the 
amplitude of the induced component is small there. As 
a result, the winding structure is reduced to the same 
configuration as that of Fig. [|(b) at higher temperature. 
Therefore, induced components have the same structure 
near T c both for P- and N-vortices, while their struc- 
tures are different at low temperature. We have con- 
firmed these winding structures also by solving the two- 
component. GL theory or the quasiclassical Eilenberger 
equationHj 

The pair potential structure for the sin (p x + p y ) ± 
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isin (— p x +p y ) case is also shown in Figs. ^] and ||. As 
for the dominant component in Fig. [2J the difference be- 
tween A' + and A'_ is smeared, as seen from the contour 
lines. The amplitude of the induced component in Fig. 

is very small, compared with the svnp x ± isinpj,-wave 
case. In this case, the gap amplitude is suppressed along 
the 0° direction on the Fermi surface. Hence, |A'_| in 
Fig. |^(c) is sharply suppressed along this direction. In 
the winding structure of the induced component for A'_ 
in Fig. ||(c), the —1- winding at the vortex center and 
the — 2-winding at the corner is the same as in Fig. [I|(a). 
But, the +l-winding points around the vortex core is ro- 
tated by 45° around the vortex center, and is located at 
the horizontal and vertical direction at low temperature. 
At higher temperatures, the +l-winding points approach 
the boundary of the unit cell, and the winding structure 
is reduced to the same as in Fig. ||(b). The transition 
temperature of the winding structure is shifted to lower 
temperatures, compared with the sinp^ ± isinp y -wave 
case. The winding of A^ in Fig. ||(d) has, however, the 
same structure as Fig. ||(b). 





FIG. 4. Winding structure of the induced component of 
the pair potential in the smp x ± isinp H -wave case, (a) The 
case of P- vortex at low temperature, (b) The case of N- vortex. 
At higher temperature, P-vortex also has this configuration. 
The plotted area is the same as in Fig. ^. The solid line shows 
the boundary of the unit cell for the square vortex lattice. The 
vortex center is presented by the solid circles. The induced 
component has the winding point at the vortex center and the 
points of the open circles. We also show the winding number 
at these winding points in the figure. In (a), total winding 
number within the dotted circle around the vortex core is +3. 
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FIG. 5. Spectrum of the LDOS N(E, r) for P-vortex 
(a) and for N-vortex (b) in the sinp x ± isinp M -wave case, 
(c) and (d) are, respectively, for A+ and A'_ in the 
sin (p x + p y ) ± isin(— p x +p H )-type pairing. The top panel 
at each figure is for the site just on the vortex center. The 
middle panel is for the site next to the vortex center. The 
bottom panel is the farthest site from the vortex center, i.e., 
the mid-point between next NN vortices. 



V. ELECTRONIC STATE AROUND VORTICES 

Since the superconducting gap of chiral p-wave su- 
perconductor opens a full gap, low energy quasiparti- 
cles are bound states around the vortex core. Therefore, 
the energy levels are discrete as in s-wave superconduc- 
tors. However, the energy levels appear at integer points. 
E n ~ nE& (n is an integer) in chiral p-wave caseJI9€3 
while levels appear M_aalf-integer points E n ~ (n+^)E& 
in the s-wave case.EZIH Here, Ea is the level spacing of 
the order Ay/Ep, which becomes narrower at higher en- 
ergy (Ao is a superconducting gap at zero field and E-p 
is the Fermi energy). In Fig. g[ we show the spectrum of 
the LDOS at selected sites inside and outside the vortex 
core. There appear some peaks within the superconduct- 
ing gap. The energy levels are found to be E\ ~ Q.\2t 
for the P-vortex in Fig. ||(a), and E\ — 0.08i for the 
N-vortex in Fig. ||(b). We observed also that the peak at 
Eq ~ vanishes just on the vortex center for P-vortex. 
Instead, sharp peaks appear there at E\ and higher E. 
The peak at E shoots up sharply at the sites next to 
the vortex center. There, the peak at E\ becomes lower, 
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and, moreover, peaks at — Ei and lower E appear. For 
N-vortex, a sharp peak is found at Eq just on the vor- 
tex center. It becomes lower at sites next to the center, 
and other peaks at iE\, ±i?2> • • ■ emerge in addition. 
Outside the vortex core, the spectrum is reduced to that 
of the zero field case. But, small peaks remain at each 
energy level in this vortex lattice case due to the low 
energy quasiparticle transfer between different vortices. 
Note also, that the broad gap edge at -E/t = 0.3 ~ 0.6 is 
due to the anisotropy of the superconducting gap. 

The energy level appearing at the vortex center is |Gep-i 
sistent with results of the single vortex calculationJi3>t3 
For an isolated vortex with circular symmetric structure, 
the energy level E e — E n can be labelled by the wind- 
ing number of the wave function u c (r) around the vor- 
tex center. For the P-vortex and N-vortex, the level E n 
corresponds to u e (r) ~ e 1 ^ 1- ™^ and u e (r) ~ e m ^, respec- 
tively. On the vortex center, u e (r) vanishes except for 
the winding state. Therefore, the LDOS at the vor- 
tex center comes from the winding state of the wave 
function. This corresponds to the level E\ (E ) for the 
P-vortex (N-vortex). While E state for the P-vortex 
is a bound state in the vortex core region, the LDOS 
vanishes at the vortex center, since the wave function is 
given by u e (r) ~ e 1 ^. In our calculation for the tight- 
binding model and the vortex lattice, the circular sym- 
metry is broken. Therefore, the wave functions with var- 
ious winding numbers are mixed for each energy level. 
Even in this case, the characteristics of the energy level 
appearing at the vortex center suggested in the single 
vortex calculation are qualitatively preserved. The mix- 
ing of the wave functions with different winding numbers 
may become larger at higher energy, because the wave 
functions extend to the boundary of the vortex lattice 
unit cell. Hence, the winding character becomes more 
smeared for higher energy state at n > 1. 

Next, we turn to the spectrum for the =p3-winding case 
of A'_j_ in the sin (p x + p y ) ± i sin (— p x + p y )-wave super- 
conductivity. It is presented in Figs. |](c) and 0(d). At 
the vortex center, the LDOS has sharp peaks at positive 
energies for the +3 winding case, shown in ||(d). For the 
—3- winding case ||(c), peaks appear at the negative ener- 
gies at the vortex center. Following the same analysis of 
Refs. |l0| and [Hj, we see that the level E n corresponds to 
the wave functions u e (r) ~ e'^ 2- ™^ and u e (r) ~ e'(' l+1 ^ 
for the +3 and —3 winding cases, respectively. Exactly 
on the vortex center, the LDOS has a peak at Ei for the 
+3-winding, and at E-\ for the — 3-winding. This wind- 
ing character affects the results of Fig. ||(c) and (d) . The 
mixing of the energy levels smears this winding character 
at higher energy levels. 

Figure || shows the temperature dependence of the 
lower five energy levels. The quantized energy levels ac- 
quire a width by forming band structure due to the inter- 
vortex transfer. Broader band width means larger trans- 
fer. The energy level at Eq ~ splits into two band for 
E > and E < 0. We show the energy levels of the up- 
per E state in Fig. [| Near T c , the levels are distributed 



with same level distances, which are small. Upon lower- 
ing of the temperature, the level distance increases. The 
level distribution then becomes different depending on 
the winding structure, as shown in Figs. ||(a)-(d). The 
increasing level distance is due to the shrinking of the 
vortex core radius on lowering temperature. In the BdG 
theory, the vortex core radius shrinks until the order of 
the atomic distance in the clean limit due to the Kramer- 
Pesch effect. Oe2i The level distance becomes larger at low 
temperature. Near zero temperature, this shrinking and 
level shift saturated. This is-,<pnsistent with results of 
the single vortex calculation. EJ'Ea 




FIG. 6. Temperature dependence of the energy lev- 
els for the P-vortex (a) and the N-vortex (b) in the 
sinpa; ± isinpy-type state. (c) and (d) are, respec- 
tively, for the and A'_ superconductivity cases for the 
sin (jp x + p y ) ± isin(— p x + p y )-type pairing. We show the 
lower five levels. Each level has energy width due to the band 
structure by the inter-vortex quasiparticle transfer. 

Since the low-energy electronic state around the vor- 
tex core determines the low-temperature behavior in the 
mixed state, we intend to further investigate the char- 
acteristics of the Eq state. The spatial structure of the 
LDOS N(E ~ 0, r) is presented in Fig. [7], where the 
contribution of i?o-states are integrated. It shows a four- 
fold symmetric star shape. Outside the vortex core, 
N(E ~ 0,r) extends toward 0° direction (i.e., a-axis and 
6-axis directions). In this direction, there is an island 
of finite LDOS around the midpoints of NNN vortices. 
In the 45° direction, the LDOS is suppressed around 
the mid-point between NN vortices. This LDOS struc- 
ture extending toward the 0° direction can be understood 
from the view point of the quasiclassical theoryclTEl In 
this framework, the low energy quasiparticle distribution 
around the vortex core is determined by the directional- 
dependent coherence length £(0k F ) noted in Sec. Ill . 
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Since £(0/c F = 0) is larger than £(0fc F = 7r /4), the quasi- 
particle propagating from the vortex core to 0° direction 
can extend farther compared with 45° direction. Thus, 
the LDOS extends toward the 0° direction. 

At the vortex center, the LDOS vanishes in Fig. f?](a) of 
the P-vortex case, reflecting the winding structure of the 
wave function u e (r) ~ e^ 1 " 71 ^ with n = 0. The LDOS 
is suppressed along the 45°-direction in the vortex core 
region around the vortex center. In Fig. ^(b) of the N- 
vortex case, the wave function it e (r) ~ e m ^ with n = 
does not vanish at the vortex center, the LDOS has a 
peak at the vortex center. It is suppressed along the 
0°-direction around the vortex center in a small region 
around the vortex core. 

The LDOS for the sm(p x + p y ) ± i sin (— p x + p y )-wave 
case is shown in Fig. CTc) and (d). Also in this case, 
the LDOS extends toward the 0° directions outside of 
the vortex core. But, the contrast between the 0° di- 
rection and the 45° direction is smeared in this case. It 
comes from the behavior of £(0fe F ). For sm(n x + p y ) ± 
isin (— p x +py)-wave case, as shown in Fig. 0(c), the ra- 
tio £(#fc F = 0)/£(6>fc F = 7r/4) is small compared with the 
smp x ± isinpy-wave case. Moreover, £(#fc F ) shows a flat 
behavior in a wider range of 8^ near minimum point 
tt/4. 

In Fig. 0(c) of the — 3-winding case (in Fig. 0(d) of the 
+3-winding case), The LDOS vanishes at the vortex cen- 
ter, reflecting the winding structure of the wave function 
u e (r) ~ e i ( n + 1 )^ (u e (r) - e^ 2 "™^) with n = 0. Since the 
wave function e 2l< ^ needs a longer distance to recover of 
the order r 2 around the singular winding point compared 
with the order r of the wave function e 1 ^, the region of 
suppression around the vortex core is wider in the +3- 
winding case. For the —3-winding case, the LDOS in the 
vortex core region is suppressed along the 45° direction. 

These LDOS structures are for the site-centered vortex 
case. For the plaquette-centered vortex, the LDOS struc- 
ture is changed within the vortex core region. It is due 
to the trapping effect of the vortex within the plaquette. 
The peak at Eq for P-vortex does not vanish near the 
vortex center, since the atomic site is not located exactly 
on the vortex center. In the both cases of the vortex 
center, we obtain qualitatively the same LDOS structure 
outside the vortex core. 




FIG. 7. Density plot of the zero-energy LDOS N(E ~ 0, r) 
for P-vortex (a) and N- vortex (b) in the sinp x ± isinp^-wave 
case, (c) and (d) are, respectively, for the A+ and A'_ super- 
conductivity cases in the sin (p x + p y ) ± i sin ( — p x + p y )-wave 
pairing. The plotted area is the same as in Fig. [| 




FIG. 8. Electron number density n(r) around the vor- 
tex at T = O.lt for P-vortex (a) and N- vortex (b) in 
the sinp x ± isinp M -wave case. (c) and (d) are, respec- 
tively, for the and A'_ superconductivity cases for 
sin (p x + p y ) ± i sin (— p x + p H )-type pairing. The plotted area 
is the same as in Fig. ^ 

Figure || shows the electron number pdansity n(r) 
around the vortex. The screening effectliS is not in- 
cluded in our calculation. The number density is 
slightly changed at the vortex center, where supercon- 
ducting.|gap-js suppressed. This is the vortex charging 
effect .E1e3 _ eZI Near T c , n(r) is slightly suppressed for all 
winding cases. With lowering temperature, the struc- 
ture is more eminent and depends on the Cooper pair 
winding. The suppression is small for N-vortex com- 
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pared with that for P-vortex. This is consistent with 
the results of the single vortex calculation.!!^ For the 
sin (p x +p y ) =p i sin (—p x + p y )- wave cases, n(r) is sup- 
pressed (enhanced) for the +3 (— 3)-winding case of the 
A'_ (Aj_) state. This n(r) structure at the vortex center 
in Fig. H is related to the asymmetry between E > and 
E < in the spectrum shown in Fig. ||. When n(r) is 
largely suppressed [(a) and (d)], the spectrum has sharp 
peaks at positive E at the vortex center. If n(r) is en- 
hanced [(c)], the spectrum has sharp peaks at negative E 
at the vortex center. When the spectrum has sharp peaks 
at E ~ at the vortex center, the suppression of n(r) is 
small [(b)]. The detailed explanation of this mechanism 
was given in Refs. |o| and 

At low temperatures, since n(r) shows the Friedel oscil- 
lation in the length order of the atomic lattice constant, 
n(r) is enhanced around the vortex along the 0° direc- 
tion before being reduced at the vortex center. For the 
sinp x =p i sin p y -wave case, n(r) exhibits a stronger sup- 
pression along the 45° direction around the vortex core, 
compared with the 0° direction. For N-vortex, since the 
suppression along the 45° direction extends farther out- 
side the vortex core, n(r) is reduced on the lines con- 
necting NN vortices. This suppression remains up to 
higher temperatures, while the suppression at the vor- 
tex core is smeared with increasing temperature. For 
P-vortex, n(r) is almost constant outside the vortex. For 
the $m(p x + P y ) ± isin(— p x +p y )-wave case, the distri- 
bution of n(r) extends toward the 0°-direction from the 
vortex core region instead of the 45°-direction. 

VI. SUMMARY AND DISCUSSIONS 

We have investigated the vortex structure in chiral p- 
wave superconductors by the BdG theory in an extended 
tight-binding model. The sinp x + isinpytype super- 
conductor and the s\np x — isinp y -type superconductor 
have different vortex structures, which are classified as 
P-vortex and N-vortex. Around the vortex core in the 
chiral superconductor, the opposite chiral component is 
induced. The induced component of the P-vortex case 
and the N-vortex case show different winding structure 
at low temperature. The winding structure of P-vortex 
depends on temperature, and it is reduced to that of the 
N-vortex case near T c . This difference in winding struc- 
ture affects on the amplitude distribution of the induced 
component and the vortex core shape of the dominant 
chiral component pair potential. In spite of the broken 
circular symmetry around the vortex core for our tight- 
binding model and the vortex lattice, we obtain qualita- 
tively the similar characteristics of the low energy bound 
states around the vortex core as thaLpf the single vortex 
calculation in a circular symmetry I 1 °H 15 I For example, we 
confirmed that the energy levels appearing at the vortex 
center and the vortex charging effect are different be- 
tween P-vortex and N-vortex. In the vortex lattice state, 



the quantized energy levels of bound states around the 
vortex core have an energy width by forming a band due 
to the inter-vortex transfer. 

We have also investigated the sin (jp x + p y ) ± 
i sin (— p x + p y ) type superconductor by considering NNN 
site pairing instead of the NN site pairing. In this case, 
the Cooper pair winding becomes =f3 at the Fermi sur- 
face. This winding structure affects the electronic state 
around the vortex core, such as the energy level appear- 
ing at the vortex center, and the charging effect at the 
vortex core region. 

In the approach of the BdG theory, we have to consider 
the case of a large superconducting gap compared to the 
Fermi energy, because of the system size is restricted in 
the numerical calculation. Further, in the clean limit, 
the vortex core radius shrinks to the order of the atomic 
distance on J^wering temperature due to the Kramer- 
Pesch effect. Thus, the level distance of the quantized 
quasiparticle around vortex core becomes larger at lower 
temperature. Since we are interested in the exotic char- 
acteristics of the quasiparticle at each quantized level, 
we calculated this case. In Sr2Ru04, this level distance 
is small. Then the contribution of the quantized level 
appears at very low temperature. We are also perform- 
ing the calculation of the quasiclassical theory, where the 
level distance is treated as zeroBj There, qualitatively the 
same structure was obtained as the result of this paper 
for the induced pair potential structure and the LDOS 
structure outside the vortex core. 

In the chiral p-wave superconductors, it is possible to 
realized a domain structure consisting of p x + ip y -wave 
and p x — rp y -wave superconducting regions. Therefore, it 
is important to clarify the difference of the vortex struc- 
ture between P-vortex and N-vortex for the investigation 
of the mixed state in Sr2Ru04. The spatial structure of 
the LDOS around the vortex include the information of 
the anisotropy in the Fermi surface structure and the su- 
perconducting gap. The observation of the LDOS struc- 
ture is important to get the information of the supercon- 
ducting gap structure and the orbital dependence of the 
Fermi surface. It may be. exaij ruaed by scanning tunneling 
microscopy experiments.EirE3~Eil 
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